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Abstract 



We consider time-dependent Gaussian wave packet solutions of the Schrodinger equation (with 
arbitrary initial central position, xq, and momentum, po, for an otherwise free-particle, but with an 
infinite wall at x = 0, so-called bouncing wave packets. We show how difference or mirror solutions 
of the form ip{x, t) — ip{—x, t) can, in this case, be normalized exactly, allowing for the evaluation of 
a number of time-dependent expectation values and other quantities in closed form. For example, 
we calculate {p'^)t explicitly which illustrates how the free-particle kinetic (and hence total) energy 
is affected by the presence of the distant boundary. We also discuss the time dependence of the 
expectation values of position, {x)t, and momentum, {p)t, and their relation to the impulsive 
force during the 'collision' with the wall. Finally, the xq,pq — > limit is shown to reduce to 
a special case of a non-standard free-particle Gaussian solution. The addition of this example 
to the literature then expands on the relatively small number of Gaussian solutions to quantum 
mechanical problems with familiar classical analogs (free particle, uniform acceleration, harmonic 
oscillator, unstable oscillator, and uniform magnetic field) available in closed form. 
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I. INTRODUCTION 



Closed-form wave packet solutions of the time-dependent Schrodinger equation are ex- 
cellent exemplary models to study the time-evolution of quantum systems for comparison 
to their classical counterparts. Because of their special properties, Gaussian solutions are 
possible for a few of the most familiar classical systems. Free-particle Gaussian wave packet 
solutions are standard fare in introductory textbooks in quantum mechanics, and this ex- 
ample goes back at least to Darwin The problem of a particle acting under a uniform 
force has similar Gaussian solutions which were first derived by Kennard m and occasion- 
ally appear in undergraduate- level presentations [3|]. Wave packet solutions for the harmonic 
oscillator are discussed in some textbooks (most often using propagator techniques) and, 
with a simple change of variables, can also be used to describe particles in unstable equilib- 
rium 5], giving rise to solutions which exhibit the expected exponential runaway behavior. 
Gaussian solutions corresponding to classical helical motion in a uniform magnetic field 
have also been constructed. In all of these cases, the special form of the Gaussian solu- 
tions allows for the explicit evaluation of time-dependent expectation values for position 
{{x)t, {x'^)t, '^Xt), momentum {{p)t, {p'^)t, ^Pt), and other quantities, for easy comparison 
to classical expectations. 

Another simple system, also with a clear classical analog, is given by an otherwise free 
particle, but restricted to the half-line by an infinite wall at x = 0, namely the problem 
defined by the ID potential 



In the classical case, a point particle would move freely, with constant kinetic (and hence 
total) energy, exhibiting an impulsive collision at the wall, resulting in a discontinuous 
change in momentum. 

Localized time-dependent solutions for this problem, dubbed bouncing wave packets, can 
be constructed in a very straightforward way from solutions of the free-particle problem. 
Andrews [3| has noted that simple difference solutions of the form ip{x,t) — tp{—x,t) not 
only satisfy the free-particle Schrodinger equation for all x values (if ip{x,t) does), but 
also accommodate the new boundary condition at the wall, namely that ip{0,t) = 0. This 
construction is very similar to image methods in electrostatics and has been used in numerical 




oo for X > 



for X < 



(1) 
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evaluations and visualizations of such systems jq|, as well as for discussions of wave packet 
propagation in the infinite square well j^, (with two infinite walls producing an infinite 
series of image wave functions.jSuch solutions have also been used to visualize many aspects 

ll| and we show in Fig. Q an example of the time-dependent 



of the collision with the wall 

t)p using just such a construction, illustrating the obvious interference effects between 
the ipi^yt) and ip{—x,t) terms, near the time of the 'bounce' with the wall. 
The study of general relationships between expectation values, such as 
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is straightforward for many of the familiarly treated quantum mechanical systems 
especially for well-behaved potential energy functions. On the other hand, for infinite well 
type potentials of the type considered here, the relationship to the classical force can be 
more subtle and so having exact or approximate closed-form solutions for quantities such as 
{x)t and {p)t to probe such relationships is very useful. For example, the force exerted on the 
walls in an infinite square well potential has been recently examined in a somewhat similar 
context, using more numerical methods, focusing on the time-dependent relationships in 
Eqn. Q as related to wave packet revivals and fractional revivals 3]- 

In this note, we will show that if one uses standard Gaussian solutions in this mirror 
or difference solution approach for the potential in Eqn. (0), one can also perform many 
(but not all) of the standard calculations including exact normalization of the wave packet 
solution, evaluation of many expectation values, and the calculation of the autocorrelation 
function, A{t), to obtain exact closed- form results, which is the main thrust of Sec. |n] 
This case then adds another example to the otherwise rather small pantheon of closed- 
form quantum mechanical Gaussian solutions to one-dimensional problems with classical 
analogs. Using these solutions, we will also be able to discuss the nature of the collision with 
the wall, deriving approximate expressions for the impulsive force exerted on the particle 
during the collision. We will also be able to explicitly evaluate the change in kinetic (and 
hence total) energy induced by the addition of the infinite barrier, making connections to 
earlier work on the effects of distant boundaries Ql on the energy spectrum of quantum 
mechanical systems; this example extends those results to an exact time- dependent wave 
packet solution and is also discussed in Sec. HH Finally, the mirror solutions discussed here 
can be connected to special, non-standard Gaussian free-particle wave packets in the special 
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limit when xo,Po 0, as shown in Sec. II I II 



II. "BOUNCING" GAUSSIAN WAVE PACKET SOLUTIONS 

We first review well-known textbook results for the standard time-dependent Gaussian 
free-particle momentum- and position-space wave packet solutions, for arbitrary initial Xq 
and po- These can be written in the forms 



'^{G){.Xjt) = ^ipo{x-xo)/h ^~iplt/2mh ^-{x-xo-pot/mf /2{ahf{l+it/to) ^ ^^-j 

^J y/7Tah{l + it /to) 

These solutions are characterized by 

{p)t = Po, Apt = (5) 

{x)t = Xo + {po/m)t = X{t), Axt = ^^/l + Wtof = ^ (6) 

where to = mha"^ is the spreading time and (3 = ah. This gives the familiar uncertainty 
principle product 

A.,.Ap, = 5yrTW - (5)(^) (T) 

for t » to- 

For Gaussian wave packet solutions of this type, the corresponding difference or mirror 
solution of Andrews {71 is written as 

7 , J AT [V^(G)(x,t) -^(G)(-x,t)] for X < 

I for X > 

where is a normalization constant. For initial free-particle wave packets, t), which 

are already normalized correctly (over all space), and which are sufficiently far apart in phase 
space {xo,po not too small), one would expect that N would be very close (exponentially 
so) to unity, a result which we will confirm by explicit calculation below. Once normalized, 
of course, any solution of the form Eqn. (jH)) will remain normalized for later times. 

In order to evaluate the various required integrals involved in the normalization of 
■?/'(G)(x, t), we make use of the fact that the difference or mirror solution is antisymmet- 
ric in X, so that |'?/'(G)(a;, t)p is automatically symmetric in x. This allows us, for example, 



to determine the normalization constant using the requirement 



N 
~2 



\i'(G){x),t)\^ dx 

|V^(G)(x,t) - V(G)(-X,t)|^C?X 
|V^(G)(x,t) - V^(G)(-X,t)|^dx 



(9) 



oo 
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and because the integration region can be extended (by symmetry) over all space, one can 
do the resulting integrals in closed form to obtain 

1 1 



where 



N 



+ 



Pop 

h 



Axo 



Po 



(10) 



(11) 



The integrals are, of course, done most simply for t = 0, but it is then easy to confirm that 
one obtains the same answer for t > 0, so that the wave function can be explicitly shown to 
remain normalized at all later times. While this expression is consistent with expectations 
for large values of xo,po, we stress that it is exact for arbitrary initial parameter values, even 
in the limit where Xo,Po ~^ 0, which is discussed in Sec. 11111 

The same method can be used to evaluate various expectation values involving even 
powers of variables, such as {x'^)t and {p)t. For example, one can show that 

{x^)t = i x^ \i)(G){.x,t)\'^ dx 



x^ \(l)(G){x,t)\^ dx 



(12) 



voty , Pi 



m 2 



where 
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and Zo is defined in Eqn. (fTT|) . The corresponding result for momentum is 



{P% = {p\ 



pI + 



(13) 



(14) 



We note the similarity in form between the expressions in Eqns. ()12j) and ()14|). so that the 
effect is the largest (not surprisingly) in the limit when zq = 0, namely when Xo,po 0. 



This special case also turns out to be a particular limit of a less-familiar free-particle solution 
which we discuss further in Sec. IIIII 

The result in Eqn. (|14p is also useful as it implies that the change in kinetic (and hence 
total) energy caused by the addition of the distant boundary is given by 

AE 2F{zo) 

E ^ l + 2ipo/3/hy' ^ ^ 

that is, an exponential suppression, which is consistent with more general arguments Q|. 

For expectation values of odd powers of x,p, the resulting integrals can be done in terms 
of error functions (erf((^)), but simple expressions for special cases of interest are perhaps 
more useful. For example, the classical prediction for the position is xciit) = —\X{t)\, 
where X{t) = (xo+Pot/fn); the time of the classical collision with the wall, tc, is determined 
by the condition X{tc) = and we can expand the integrals required for an evaluation of 
{x)t in terms of X{t ~ tc), treated as a small parameter, to obtain the approximation 

,2._/ A [x{t)r 



{x)t^t.= I x\^^G){x,t)\'dx^ (-^-^^^ + --- ] . (16) 



The first term in this expansion was derived in Ref. We show, in Fig. ^ the numeri- 
cally calculated value of {x)t versus t for a sample solution and note how the expression in 
Eqn. ()16|) reproduces the softened 'parabolic' shape of the curve near the collision time, to 
be compared to the purely classical result, namely xcbif) = —\X{t)\, shown as the dashed 
lines. 

The result in Eqn. fll6|) can then be differentiated and evaluated at the collision time 1^ 
to give 



d{x)t 



dt 



h 



t^tc 



(3^ 



' ' (17) 



Vi + (tcAo)' 

for tc » to- This implies that the expectation value of momentum at the classical collision 
time does not vanish, for reasons which are most easily visualized using the behavior of the 
time-dependent momentum-space probability densities, as in Fig. 3 of Ref. . For times 
long before and long after the 'bounce', the momentum distribution is peaked at +po (within 
±Apo) and —po (within ±Apo) respectively. Roughly speaking, at the classical collision time, 
the momentum components corresponding to the fastest speeds (in the +pq + Apo half of 
the distribution) have already been reflected from the wall and are already flipped in sign to 
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have values — (po + ^Po), while the slower components are still predominantly in the lower 
+Po — Apo half of the distribution, but still positive. There is a small resulting asymmetry 
in the momentum-distribution, giving an expectation value of order —ApQ, as in Eqn. (|17|). 

The result in Eqn. is also useful in that it gives similarly valid approximations for 
higher derivatives, such as 



m- 



(18) 



which can then be used to describe the effective force exerted on the particle by the wall near 
the classical collision time. It is not clear how one would obtain this result more directly from 
the Ehrenfest theorem approach in Eqn. (0) using the potential energy function in Eqn. (^J. 
The dimensional dependence of this result on ^) can be easily understood from simply 



assuming that the change in momentum during the collision is of order Ap = Pf—pi 
while the collision time. At, is determined by 



Po 



m 



Ax 
At 



or 



At 



givmg 



Ap 
At 



-2po, 



(19) 



assuming we identify Ax ~ (3t near the time of the collision 

Finally, another useful quantity which can be evaluated in closed form for the solution 
in Eqn. (jSJ is the autocorrelation function, which measures the overlap between the initial 
quantum state, '?/'(x, 0), and its time-evolved value at later times, ipi^yt)- This is defined 
most generally by 



A{t) 



{x,0)]*ip{x,t) dx 



For the free-particle Gaussian solutions in Eqns. Q and 
equivalent forms [3| 



[0(p,O)]*0(p,t)rfp. (20) 
) this can be written in the 



^(G)(t) 



^1 - zt/2t( 
1 

a/1 - zt/2tt 



{x{t)-xoY 

4/^2(1 -2t/2to) 



; exp 



ia^plt 



2to(l - it/2to) 



(21) 
(22) 



Both of these forms give 



+ (V2to) 



: exp 



~2a'pl- 



{t/2tof 



l + (t/2to)2) 



(23) 
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One sees that the free-particle autocorrelation function decreases monotonically with time, 
due to both the dynamic exponential dependence on po? as well as to the dispersive pre-factor 
which can be attributed to wave packet spreading. 

For the bouncing wavepacket solution, we must evaluate the autocorrelation using the 
closed form expression for the position-space solution in Eqn. (jS)), giving 

^(G)it) = / [i'(G){x,0)]*tjjG){x,t)dx 



oo 





= / [^(G)(x,0)-7/;(G)(-x,0)]*[V'(G)(a;,t)-V'(G)(-a:,t)]cia; (24) 

J —oo 

and the same trick of extending the integral over all space (because the integrand is still 
odd under x —x) can be used to find 

■1 - exp [- {{xo/PY + ipoP/hf} /(I + it/2to)]' 



AiG){t) = A(G)(t) 



1-exp [-xo/f3y + {pof3/hy] 



- A,a,it) ^[-°/(^ + f^°)n . (25) 

Once again, there is a monotonic decrease in \ A{t) \ with no distinction between the smoother 
time-evolution before and after the collision and the time during the impulsive splash at the 
wall. 



III. OTHER FREE-PARTICLE GAUSSIAN SOLUTIONS RELATED TO THE 
BOUNCING WAVE PACKET 

While the standard Gaussian solutions for the free-particle case given in Eqns. and 
(0} are the most familiar examples found in textbooks, it is straightforward to construct 
other localized Gaussian-like wave packet solutions, using the fact that a wide variety of 
Gaussian integrals can be performed in closed form. Some of these can then be easily used 
as special case solutions for the 'bouncing' wavepacket well. 

For example, a modified free-particle momentum-space solution of the form 



gives the expectation values 



{p)t = Po, {p^)t = pI + and = ^ ^ (^7) 



and can be Fourier transformed to yield the position-space wavefunction 



t) = 2 J ^3^^ ^ ^^^^ ^3 e"P'^(--o)/n^-^plt/2mn _ ^^^^^ ^_(,._x(f))V2/3^(l+.tAo) (28) 



with 



{x)t = X{t)=xo + ^— and Axi = w|A. (29) 



The uncertainty principle product in this case is given by 

/\x,-Ap, = —^lTWJW (30) 

which is similar to that for the first excited state of the harmonic oscillator, at least for 
t = 0. In the same way, initial momentum-distributions with higher powers of (p — Po) can 
be used to exhibit localized position-space wave packets. 

For the case described by Eqns. (j26|) and (j28|l . if we consider the special case of xo,po 0, 
we find the solution 



which is valid for all space. This clearly satisfies '?/'(0,t) = for all t and so can be used 
as a solution for the 'bouncing' packet case corresponding to the potential in Eqn. (P); the 
solution must then be 'renormalized' by multiplying by a factor of \/2 to account for the 
different range of definition. This then gives a 'bouncing' packet solution 

. (32) 

for a; > 
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The 'bouncing' wave packet solution of Eqn. (fTn|) which is valid for arbitrary values of 
Xq and Po can be considered in the limit when = and po ~^ 0. In that limit, we recover 
the form in Eqn. ()32|) (except for trivial multiplicative factors of i.) 

For this special solution, almost all of the relevant expectation values can be obtained in 
closed form, and we find, for example, that 



(x), = -^ and (x^), = ^ (33) 



which combine to give 



A / Stt-S 
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which increases with time in the same way as for the standard free-particle solution in 
Eqn. ®. 

Using the explicit form for ip(_Q-^{x,t) in Eqn. (j32p . and the operator form for p = 
{h/i){d/dx), we also find that 



2h 

{p)t = 



t/to 



and = 1^ . (35) 




(t/to) 

The first of these expressions is consistent with differentiation of the result in Eqn. while 
the second is related to the free-particle result in Eqn. (jTfj) . The momentum uncertainty is 
then given by 

+ it/tor] ^ ^ 

which actually decreases in time. This effect can be understood crudely as being due to 
the fact that the t = solution has both positive and negative momentum components 
in the range {—Apo,+Apo), and over a time interval of order to, the positive momentum 
components are reflected from the inflnite wall so that the momentum distribution is then 
more localized to the range {—po,0). 

The effective force on the particle due to its interaction with the wall, can be associated 
with 

^ = - (^-] ^ (37) 

which decreases monotonically with time, with an initial value which scales as Apo/^o which 
is dimensionally correct. 

These results can be combined to give the uncertainty principle product as 



h /3(37r-8) 



Axt ■ Apt = - W ^ ^ Vl + (1 - 8/37r)(t/to)2 (38) 



(0.58n)v/l + (0.15) (t/to) 



2 



This localized solution has an initial uncertainty principle product which is only slightly 
larger than the minimum value but for long times is substantially smaller (about half) 
than that of the standard Gaussian free-particle solution in Eqn. ((7j), since for t/to >> 1; 
Eqn. (jHHj) reduces to 

A,..Ap, _ ^ . (0,45)1 fiV (39) 



2\ TT J \to J 2\to 
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Variations on the momentum-space distribution in Eqn. ()26|) which are also odd in p — po 
(higher power solutions of the simple harmonic oscillator, for example) can be used to 
evaluate other generalized solutions for the free-particle case which can also satisfy the 
appropriate boundary conditions (at a; = 0) for the 'bouncing' particle case and provide 
additional examples. 

IV. CONCLUSIONS AND DISCUSSIONS 

Using the familiar method-of-images or mirror wavefunction technique, we have shown 
how to construct normalizable Gaussian solutions with arbitrary initial xq,po for the bounc- 
ing particle problem, for which many expectation values and related quantities are calculable 
in closed form. This example adds another case to the limited number of time- dependent 
wave packet solutions of one-dimensional quantum mechanical problems with familiar classi- 
cal analogs. It also provides an explicitly calculable example of the effect of the introduction 
of a distant boundary on an explicitly time- dependent solution for a quantum system for 
comparison to more general discussions. Because of the methods used, not all of the ex- 
pectation values (or related quantities such as the momentum-space wave function or the 
Wigner quasi-probability distribution) can be evaluated as easily, but all of the even ex- 
pectation values are readily calculable and many of the others, such as {x)t can be usefully 
approximated, especially near the time of the classical collision with the wall. A class of 
non-standard free-particle Gaussian solutions can also be used to provide special pq = 
which vanish at x = and are therefore also solutions for the bouncing well case, with 
interesting long-time behavior. 
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l\|/(x,t)P vs. X t vs. <X>j 

FIG. 1: Plots of the position-space probability density, \^p(^Q•^{x,t)\'^ versus x, for the Gaussian 
bouncing packet of Eqn. 0, for various times before, during, and after the 'collision' with the wall 
are shown on the left. On the right, we show numerical calculations of the quantum mechanical 
expectation value of position, {x)t versus t (solid curve), over the same time range. The softening 
of the classical trajectory result, xcL^t) = —\X{t)\ (dashed lines), near the collision in the quantum 
case is well described by Eqn. (fT6|) . 
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